Nonmonotonic energy harvesting efficiency in biased exciton chains 
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We theoretically study the efficiency of energy harvesting in linear exciton chains with an energy 
bias, where the initial excitation is taking place at the high-energy end of the chain and the energy 
is harvested (trapped) at the other end. The efficiency is characterized by means of the average 
time for the exciton to be trapped after the initial excitation. The exciton transport is treated as 
the intraband energy relaxation over the states obtained by numerically diagonalizing the Frenkel 
Hamiltonian that corresponds to the biased chain. The relevant intraband scattering rates are 
obtained from a linear exciton-phonon interaction. Numerical solution of the Pauli master equation 
that describes the relaxation and trapping processes, reveals a complicated interplay of factors 
that determine the overall harvesting efficiency. Specifically, if the trapping step is slower than or 
comparable to the intraband relaxation, this efficiency shows a nonmonotonic dependence on the 
bias: it first increases when introducing a bias, reaches a maximum at an optimal bias value, and 
then decreases again because of dynamic (Bloch) localization of the exciton states. Effects of on-site 
(diagonal) disorder, leading to Anderson localization, are addressed as well. 
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I. INTRODUCTION 

The harvesting of electromagnetic energy, i.e., its ab- 
sorption, transport to a specific site, and subsequent 
trapping in an alternative form of energy, is a process of 
great importance to life on earth. In nature, systems oc- 
cur that perform these tasks with amazing efficiency. In 
order to make the entire process unidirectional, a certain 
downward energy gradient (bias) occurs. An excellent ex- 
ample is the light-harvesting system in purple bacteria, 
where photons with a wavelength around 800 nm may be 
absorbed in the B800 ring of the light-harvesting system 
2 (LH2). The resulting excitation may be transferred to 
the lower-energy B850 ring (850 nm), before undergoing 
transfer to the yet lower-energy LH1 system (875 nm), 
from where the final transfer to the so-called special pair 
of the reaction center occurs ("the trap"). Here, the ex- 
citation energy is transformed into chemical energy (see 
Refs.[l|andUfor a review). All these steps occur rapidly 
and with extremely high efficiency. It is well-understood 
by now that the transport in this system is not described 
by the classical Forster mechanism^ of incoherent en- 
ergy transfer. The excitonic nature plays an important 
role, both for transport inside and between the various 
subsystems ,2i5jSiZiS 

Also in man-made nanoscale systems, the injection of 
energy at some point and its subsequent transport to and 
collection at another site, is a process of growing inter- 
est. Not only is this relevant to artificial light-harvesting 
systems, but also to electronic and photonic nanodevices. 
Obviously, it is important to study how the efficiency of 
the total process can be optimized. 

The aim of the present paper is to study theoretically 
the effect of a downward energy gradient (bias) on the ef- 
ficiency of energy harvesting in linear chains of strongly 
coupled two-level units, and thus to investigate how this 
process may be optimized. In the context of this paper, 



we consider harvesting as the total process following ini- 
tial excitation at the high-energy side of the chain until 
its trapping at the low-energy side, i.e., it includes both 
the transport and the trapping processes. We will choose 
the average trapping time as a measure for the efficiency 
of the harvesting process. 

At first glance, it seems obvious that a downward bias 
will reduce the harvesting timescale, and thus will in- 
crease its efficiency. Indeed, in the Forster limit of in- 
coherent hopping between weakly coupled sites, an en- 
ergy bias will cause the excitation to diffuse towards the 
lower energy end of the chain and will generally lead to a 
higher quantum efficiency for the harvesting process. Ex- 
perimental realizations of such chains have been studied 
by Sauer and coworkers j^i in addition low-generation 
dendritic systems with an energy bias towards their core 
have been studied by various groupsj 11 ' 12 ! 13 

In contrast to the weak-coupling case, the effect of a 
bias is not clear a priori in chains with strong intersite 
coupling, where in the absence of a bias the excited states 
are delocalized (excitons) over the entire chain. On the 
one hand, the bias defines a preference for the excitation 
to move towards the trap, on the other hand, it leads to 
dynamic (Bloch) localization of the exciton stateS ) 14 i 15 
which slows down the transport process. We will show 
that, as a result of the competition between both effects, 
a strong bias always decelerates the harvesting process, 
while for smaller values of the bias, nonmonotonic behav- 
ior may arise. We also study the effect of energetic disor- 
der, which causes Anderson localizatio n 16 ! 17 and thereby 
also acts toward reducing the transport efficiency. 

To model the harvesting process, we will describe the 
system by a Frenkel exciton Hamiltonian, with site ener- 
gies that follow a linear energy bias and, superimposed 
on that, may suffer from random disorder. The initial 
state is a localized excitation of a single monomer at the 
high-energy end of the chain, while at each time the rate 
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of trapping is proportional to the probability that the 
monomer at the low-energy end of the chain (the trap) is 
excited. Thus, to get trapped the exciton should either 
overlap appreciably with the trap or travel over the chain 
until it does so. The transport is modeled as an intra- 
band relaxation process, induced by the scattering of the 
exciton eigenstates on phonons, and described by a Pauli 
master equation for the exciton populations. We will re- 
strict ourselves to the zero-temperature limit. A study 
of temperature effects will be presented elsewhere^ 

Possible realizations of our model are molecular ag- 
gregate a 19 ' 20 and conjugated polymer a 21 i 22 subjected to 
a non-uniform electric field in order to create a bias 
(Stark effect), molecular photonic wires 9 - ' 10 ' 23 in which 
a bias is created by chemical synthesis of a series of 
molecules with different transition energies, linear arrays 
of metal nanoparticles , 24 ' 25 and assemblies of semicon- 
ductor quantum dots i 26 ' 27 ' 28 In the latter two examples, 
a bias may be created by carefully arranging particles of 
different size. Another possible application of the under- 
lying physics might be the electrical transport through a 
single DNA molecule j 29 ' 30 ' 31 ' 32 the source-drain voltage 
applied to the DNA is usually on the order of a few volts, 
which is high enough for dynamic localization to occur. 

The outline of this paper is as follows. In the next 
section, we present our model Hamiltonian. The Pauli 
master equation for the population dynamics as well as 
the trapping model are introduced in Sec. lIIII In Sec. lIVi 
we present our numerical results, analyzing first the effect 
of a bias on the exciton trapping time in disorder-free 
chains (Sec. IIV Al) . The interplay of diagonal disorder 
and a bias is unraveled in Sec. IIV Bl We summarize and 
make final remarks in Sec. |Vj 
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FIG. 1: Schematics of the system under consideration: a one- 
dimensional array of two-level units with a linearly varying 
excitation energy. Static disorder can be present as well, in- 
troducing an additional stochastic variation in the excitation 
energies. The initial excitation is created at the high-energy 
site n — N. The excited state of site n = 1 is strongly 
quenched by an additional molecule, so that it acts as a trap. 



II. MODEL HAMILTONIAN 

We consider a one-dimensional chain of N two-level 
units with aligned transition dipoles, which are coupled 
through dipole-dipole transfer interactions. The site ex- 
citation energies, in addition to a stochastic component 



(diagonal disorder), include a linearly varying part (bias). 
The schematics of the system is depicted in Fig. [TJ The 
low-energy end of the chain contains a trap, while en- 
ergy is injected at the high-energy end. In Sec. IIII1 we 
will define how the transport and trapping processes take 
place. In this section, we consider the relevant Hamilto- 
nian, which reads 

N N 

H ex = 5^[e„+(n-l)A] \n) (n|+ J n m\n)(m\, (1) 

n—l n,m^£n 

Here, the state \n) corresponds to site n being excited, 
while all other sites are in the ground state. The site 
excitation energy contains two parts: (i) £„, which is 
a Gaussian stochastic variable with mean e (which can 
be set to zero without loss of generality) and standard 
deviation a, and (ii) the energy bias, characterized by 
the parameter A. The transfer integrals J nrn are given by 
Jnm = ~J/\ n ~ m \ 3 i with the nearest-neighbor coupling 
-J < 0. 

The exciton eigenenergies E s (s = 1, ...,N) and eigen- 

functions |s) = J2 n =i Csn I")' with c sn being the nth com- 
ponent of the sth eigenfunction, are found by diagonaliz- 
ing the matrix H nm = (n\ H \m) = 5 nm [e n + (n- 1)A] + 
J nm and can be taken real in our model. In the disorder- 
and bias- free case (a = A = 0), the exciton states are 
delocalized over the whole chain. As a convenient refer- 
ence case, we note that the wave functions |s) resemble 
those obtained when only the nearest-neighbor coupling 
is taken into account! 33 ' 34 ' 35 

/ , \ 1/2 IV / \ 

x ' n—l N ' 

The N exciton states form a band with approximate en- 
ergies 

N i / \ 

n—l v 7 

that range from E x = -2.404J to E N = 1.803J, thus 
having a width E N - E x = 4.207 JM Note that E s in- 
creases with s (J > 0). Throughout this paper (i.e., also 
in the presence of bias and/or disorder), we will use the 
convention that s labels the states in order of growing en- 
ergy. Although Eqs. and © are useful for reference, 
it should be stressed that all results in this paper are 
obtained by exact numerical diagonalization, accounting 
for all dipole-dipole interactions. 

In a disorder-free system, a nonzero bias results (in the 
thermodynamic limit N — > oo) in dynamic (Bloch) local- 
ization of all states! 14 ' 15 The localization size is estimated 
from semiclassical arguments as Lb = {Em — Ei)/A. 36 
The Bloch localization is accompanied by the subsequent 
reorganization of the energy spectrum of the system, 
which becomes ladder-like with level spacing Aj 3 - 7 - This 
structure is revealed in photoluminescenc o 38 ' 39 and pho- 
toconductivity^ - spectra of semiconductor superlattices 
as a series of equally spaced peaks. 
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FIG. 2: Exciton wave functions of a disorder-free chain of 
size N — 50 for the states s = 1 (dashed curves), s = 5 
(dash-dotted curves), and s = 25 (solid curves) calculated by 
diagonalizing the Hamiltonian {1} for various magnitudes of 
the bias: A = (upper panel), A = 0.1J (middle panel), and 
A = 0.3 J (lower panel). The lines are a guide for the eye; the 
coefficients c s „ are only defined for integer values of n. 

In finite-size systems, dynamic localization comes into 
play when Lb becomes smaller than the system size N 
or, in other words, when the overall bias NA exceeds 
the bare exciton bandwidth 4.207 J. Figure [21 where we 
plotted several wave functions calculated for three mag- 
nitudes of the bias, illustrates this. Note that for Lb < N 
there exists a strong correlation between the energy of an 
exciton and its location: the states with smaller s (near 
the bottom of the exciton band) are localized at the lower 
energy end of the chain, while the states with higher s 
(near the top of the band) are localized at the higher en- 
ergy end. Upon increasing A, the bandwidth approaches 
a value TV A ii 

Introducing a stochastic component in the site energies 
e n will result in an additional (Anderson) localization 
effect. At sufficiently high disorder strength a, the bias- 
induced correlations between the energy of a state and 
its location on the chain will be destroyed. The lowest- 
energy exciton states will not necessarily be located at 
the low-n side of the chain. 



III. TRANSPORT AND TRAPPING MODELS 

A. Formalism 

The proper description of the dynamics of quantum 
particles in a system that interacts with a bath, requires 
the density matrix formalismi 42 i 43 ' 44 This holds in partic- 
ular for the transport and trapping processes underlying 
the energy harvesting considered here. The system stud- 
ied is described by the reduced density matrix p ss i of the 
exciton system, which is obtained from the total den- 
sity operator by tracing out the bath. In our case, the 



bath consists of the vibrations in the host surrounding 
the exciton chain. The diagonal elements p ss represent 
the exciton populations, i.e., the probability that exciton 
eigenstate s is excited; the p ss > for s ^ s' describe the 
coherences between different exciton states. 

For sufficiently weak coupling between the excitons and 
the host vibrations, we can use the Born-Markov approx- 
imation, where it is assumed that the interaction with the 
host can be treated as a perturbation and that the host 
is always in equilibrium. A further simplification occurs 
if we only consider time scales that are long compared 
to the typical oscillation times of the density matrix ele- 
ments. This is referred to as the secular approximation, 43 
and leads to a decoupling of the time evolution of the 
populations and the coherences. 

The time evolution of the populations is then given 
by the Pauli master equation and does not involve any 
coherences: 

Pss = ^ PrrWsr ~ Pss ^ Wrs ' ( 4 ) 

r^s r^s 

where the vibration-induced scattering rates W ss i will be 
discussed later on. In the secular approximation, the co- 
herences p ss i decay exponentially with the dephasing rate 
Iss' = \ Yl<r(Wrs + W rs i) while oscillating in time at a 
rate lo ss i — E s — E s i . Additional terms will be included 
to account for the trap and spontaneous decay and, in 
fact, the trapping term does introduce a new coupling of 
the populations to the coherences, as will be seen shortly. 
However, in most cases we consider, the dephasing of the 
coherences is sufficiently rapid compared to the harvest- 
ing time to allow for a treatment that involves the pop- 
ulations only. A further discussion on the validity of this 
approach is presented in Sec. IIII Dl 

As we already mentioned in the Introduction, the en- 
ergy harvesting process not only includes the transport, 
but also the trapping at the low-energy end of the chain. 
We will assume that the trapping does not affect the ex- 
citon wave functions and that this process changes both 
the exciton populations and coherences with a rate that 
is proportional to the probability that site n = 1 is ex- 

cited: p { Z? P) = -hnSm + Snn)pt? P) - We will refer to 
the constant T as the quenching or trapping amplitude. 

By transforming p^n?^ to the exciton basis, it is found 
that indeed both the exciton populations and coherences 
are affected by the trap and that the action of the trap 
couples populations to coherences. If on the timescale of 
the harvesting the coherences may indeed be neglected, 
the population of any exciton state s decays nonradia- 
tively with the rate 

r s - Tc 2 sl . (5) 

After accounting for spontaneous emission by introducing 
a radiative decay channel with rate 7 S = 7oEn c sn) 2 
(70 is the radiative decay rate of an isolated monomer) , 
the overall exciton population dynamics can be described 
by the following Pauli master equation for the exciton 
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populations P s = p ss : 

p s = -hs + r s )p s + J2 (Wss>p s > - u.j '. i . (6) 

s' 

As the initial condition to Eq. ((6|), we will assume 
that the rightmost site n — N is excited at t = 0: 
Pnn'(0) = S nn 'S n N- A transformation to the exciton ba- 
sis yields p ss '(fi) — c s nc s >n- This corresponds to a nor- 
malized initial exciton population distribution given by 
P s (0) = c 2 sN . Furthermore, our initial condition implies 
that off-diagonal density matrix elements p ss ' (exciton 
coherences) are also created; the effect of these coher- 
ences on the energy transport may be relevant when their 
dephasing times are larger than or comparable in magni- 
tude to the typical harvesting times. This is the case for 
a sufficiently strong quenching amplitude T in combina- 
tion with strongly delocalized exciton states. We come 
back to this point in Sec. IIIIDI 



B. Exciton scattering rates 

In the Born-Markov approximation, the scattering 
rates W ss ' which feature in Eq. (QJ are equivalent to those 
obtained using the Fermi golden rule, with the exciton- 
vibration interaction serving as the perturbation. As 
stated above, we think of vibrations in the chain's host 
medium as the ones mainly responsible for the scatter- 
ing. We will restrict ourselves to one-phonon processes 
and use a scattering rate from state s' to state s of the 
formed! 

N 

W ss >= Wo S{\E S -E S '\) Y,<?sn<?s>n 

n=l 




n[E s -E s ,), E S >E S ,, 
l + n(E s ,-E s ), E S <E S , 



Here, the prefactor Wo is a measure of exciton scattering. 
The spectral density S(\E S — E s > |) depends on the details 
of the exciton- vibration coupling as well as on the vibra- 
tional density of states. In particular, within the Debye 
model for the vibrational modes, this factor is given by 
S(E S -E B ,) = (\E S - E s , |/ J) 3 (see, e.g., RefEzj). We use 
a Debye-like spectral density with an exponential cut-off 
factor, S(E S -E S ,) = \(E S -E S ,)/J\ 3 exp(-\E s -E s ,\/uj c ), 
where uj c is a cut-off frequency. Similar spectral densities 
have been successfully used to fit the optical dynamics in 
photosynthetic antenna complexes. 44 ! 48 ! 49 ' 50 Note, how- 
ever, that the results that we will present later on only 
weakly depend on the particular form of S(E S — E s i). 

The term Y^n=i c L c s'ti m Eq. ® represents the over- 
lap integral of exciton probabilities for the states s and 
s' . It depends on both the magnitude of the bias, A, and 
the strength of the disorder, a. Because of the overlap 
integral, only hops between neighboring localized states 



are efficient. This strongly limits the exciton transport 
along the chain. 

Finally, n(E) — [exp(E/T) — 1] _1 is the occupation 
number of the vibrational mode with energy E (the 
Boltzmann constant is set to unity) . Due to the presence 
of the factors n(E) and 1 + n(E), the rates W ss i meet 
the principle of detailed balance: W ss i = W s ' s exp[(E s — 
E s t)jT\, Thus, in the absence of decay channels, the 
eventual exciton distribution is the Boltzmann equilib- 
rium distribution. 



C. Harvesting time 

The object of our study will be the harvesting time r 
which we define as follows. We first introduce the overall 
survival time of an exciton as the standard expectation 
value: 

f= [°° dtt(-Y J ps))= r dt(j2Ps(t)) , (7) 

Jo x ' Jo x s 1 

where the angular brackets denote the average over dis- 
order realizations. This quantity should be calculated 
for two situations: with (f) and without (jo) the trap. 
Finally, the exciton harvesting time t is obtained by ex- 
tracting the survival time to with respect to the radiative 
decay from the overall survival time f , according to the 
rule)£i 



In fact, calculating r does not require a full solution of 
the kinetic equation ([6]). Indeed, its solution reads 

A(t) = E(^)^'(0), (9) 

s' 

where we have introduced a matrix R with elements 

Rss> = (is + r s + w s „ s )s ss , - w ss , . (io) 

s" 

Using Eq. ([9]), the integration over time in Eq. ([7]) can 
be performed formally, yielding 

Thus, to calculate r for a particular choice of parameters, 
one has to invert two N x N matrices R, namely one for 
the chain with trap and one excluding the trap. 

D. Validity of the formalism 

Throughout most of this paper, only the exciton popu- 
lations are considered, which in most cases is an excellent 
approximation. The limit T ^> Wo, however, highlights 
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some limitations of our model, which are caused by not 
taking the coherences p ss > into account. In most situa- 
tions we consider, this is justified, as the coherences p ss > 
decay with a dephasing rate of an order of magnitude 
of ^2 r (W ra + W rs i). This implies that as soon as some 
relaxation process is involved in reaching the trap, the 
harvesting process will be slow enough to allow for the 
neglect of coherences. However, when the quenching am- 
plitude T is large and the exciton states are delocalized 
over the entire chain (i.e., short chains, small bias and 
small disorder), our model also produces a clearly un- 
physical artifact, namely the possibility of instantaneous 
energy transfer from the site n — N to the trap. This 
can be clearly seen from the quenching term in the Pauli 
master equation —T S P S is non-zero at t = for the 
initially populated exciton states. This problem is reme- 
died by including coherences in the trapping term, which 
in turn leads to coupling between different coherences 
and populations. 

We can estimate when the coherences may safely be 
neglected by requiring that the probability of the initial 
population distribution undergoing relaxation should be 
larger than the probability of being directly caught by 
the trap: £ s (£✓ W s , s c? sN ) > E.IV&c^. For the 
case of no bias and no disorder, we can use the (approx- 
imate) explicit expressions for c sn in Eq. ^ to numer- 
ically estimate the above inequality. The left hand side, 
which depends on the form of the scattering rate, nu- 
merically evaluates to a iV-independent value of approxi- 
mately Wo/20, while the right hand side can be estimated 
as ^jf. This estimate shows that our formalism is valid 
for T < AWo/30, which for the chain length N = 50 
and the scattering amplitude W = 10J used later on 
translates into V < 17 J. Note that at larger biases Bloch 
localization strengthens the validity of our model, as the 
delocalized exciton states extending from the initially ex- 
cited monomer to the trap that are required for direct 
quenching (i.e., with no relaxation involved) become in- 
creasingly rare and relaxation will thus be necessary. In 
that case, the coherences will dephase sufficiently fast to 
become irrelevant on the total timescale of the harvest- 
ing. 

In order to get a more quantitative impression of the 
effect of coherences, we have performed some additional 
calculations in a crude model where, within the secular 
limit, coherent effects on the trapping were taken into 
account by assuming a simple exponential decay of the 
coherences, while both population and coherence trap- 
ping terms were included: 

■(trap) = -Tc 2 slP t aP) ~ \VC SI £ C sn (p„/ + Ps's) , 

r^s 

(12a) 

Ass' = - iu} ss <p ss > - - ^2(W rs + W rs ')pss' 

r 

- \y (c 2 s1 + c 2 s/1 ) Pss , . (12b) 



These equations are obtained by disregarding all terms 
that couple the time evolution of the coherences to other 
density matrix elements, while the equation for the time 
evolution of the populations is complete and naturally 
also includes the phonon scattering and radiative decay 
terms discussed previously. Although this is not a con- 
sistent approach, it does allow for a first impression of 
the effect of coherences on population decay. As we will 
see in Sec. IIV1 the results of this approach confirm the 
above estimates. 
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FIG. 3: Harvesting time r [denned in Eq. ([H]l] versus the 
applied bias magnitude A calculated for a disorder-free chain 
of 50 sites for various quenching amplitudes T (shown in the 
plot). The exciton scattering amplitude was set to Wo = 10 J. 



IV. NUMERICAL RESULTS AND DISCUSSION 

All results reported in this paper were obtained by us- 
ing physical parameters (J, 70, and Wo) typical for a par- 
ticular physical system where our theory might be appli- 
cable, namely molecular aggregates. More specifically, 
we set J = 600 cm" 1 and 70 = 3.6 x lOV 1 = 2 x 10~ 5 J, 
which are appropriate values for J-aggregates of the dye 
pseudoisocyanine (PIC) i 52 i 53 The scattering amplitude 
Wo was fixed at Wo = 10J, which is consistent with 
the values found from the fit of the temperature depen- 
dence of the J-bandwidth and of the radiative lifetime 
in J-aggregates of PICi^ The cut-off frequency uj c of the 
spectral density S(E) was taken equal to lo c = 0.5 J. The 
exact value is arbitrary to a large extent and hardly af- 
fects the conclusions derived from our study, although it 
is crucial that some cut-off value is introduced to sup- 
press unphysically energetic phonon modes. To simplify 
the link to other systems, the nearest-neighbor coupling J 
was used as the unit for all relevant quantities. We stress 
once again that all dipole-dipole couplings are taken into 
account. Chains of size N = 50 were considered in all 
simulations. As was already mentioned in the Introduc- 
tion, we only consider the zero-temperature limit. 
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A. Disorder-free biased chain 

Figure [3J shows the behavior of the harvesting time r 
as a function of the bias magnitude A, calculated for a 
set of quenching amplitudes T, indicated in the plot. As 
is seen, the character of the behavior is strongly affected 
by r, revealing an acceleration of the harvesting process 
with increasing bias for T < Wo(— 10 J), while the A- 
dependence at T ~ Wo is shown to be nonmonotonic. 
For r 3> Wo, the model predicts a harvesting time that 
monotonically increases with A, starting from r« 1/J 
for A = 0. We do not show and analyze this limit how- 
ever, as it describes an unphysically strong quenching ef- 
fect, with r being much larger than the nearest-neighbor 
coupling J. Such strong quenching would affect the exci- 
ton wave functions; moreover, the neglect of coherences 
is not a valid approximation in this case (cf. Sec. IIII D)) . 

To unravel the behavior found in Fig. [3J we note that 
various bias-dependent factors exist that govern the over- 
all energy harvesting process: (i) - the initial distribu- 
tion of population -P s (0) = c 2 N , (ii) - the trapping of 
the exciton states with rate T s — Tc 2 sl , and (iii) - the 
intra-band relaxation which populates (or depopulates) 
the most strongly trapped states. The interplay of all 
of these factors determines the bias dependence of the 
exciton harvesting presented in Fig. [3J 

We first discuss the initial condition, P s (0) = c 2 N , 
which in Fig. 0Ja) is plotted as a function of the state 
index s for four choices of the bias strength A. In the 
absence of a bias (A ~ 0), the states in the center of 
the band are populated to a much larger extent than the 
states near the bottom and the top of it. Similar behav- 
ior occurs in the nearest-neighbor coupling model, where 
this can be explicitly shown by replacing c s jv with the 
iV-th coefficient in Eq. (J2]). A nonzero bias drastically 
changes this situation, because of the dynamic localiza- 
tion. The localization dictates that the states which are 
located near the right-most site n — N, i.e. those which 
have a significant amplitude at the site of initial excita- 
tion, are those with higher energies or large s (see Fig. [5]). 
Thus, upon increasing the bias strength, the initial dis- 
tribution of population will shift to the top of the exciton 
band. Figure Ufa) clearly illustrates this. 

The distribution c 2 sl , which determines the quench- 
ing rate T s , shows the opposite tendency, as is seen in 
Fig. 0Jb). In the absence of a bias, it is maximal for 
the central band states and negligible at the band edges, 
similar to c 2 sN . Already at a moderate bias, the quench- 
ing rate of the states at the bottom of the band increases 
significantly. For larger bias magnitudes, the distribution 
c 2 sl strongly shifts to the lower band edge (to small s). 
For a sufficiently strong bias (A A > E N -E 1 = 4.207 J), 
when the Bloch localization sets in, the initially excited 
states and those undergoing efficient quenching are sepa- 
rated completely and located at the top and the bottom 
of the band, respectively. 

Finally, an important stage in the harvesting process 
is the relaxation of excitons from the initially populated 
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FIG. 4: Bias dependence of the exciton probabilities c? aN (a) 
and (b) that determine the initial distribution of the exci- 
ton population P s (0) and the quenching rate T a , respectively. 
The bias tends to shift the initial distribution of population 
to the top of the band (to higher s). Oppositely, the strongly 
quenched states are shifted to the lower band edge (small s) . 

states to those where the quenching rates T s = Vc 2 sl be- 
come comparable to or larger than the relaxation rate 
XL' Ws's- This relaxation process will generally take 
longer when increasing the bias magnitude A. The reason 
for this is twofold. First, as we have seen above, for larger 
bias the excitation and quenching occur at opposite band 
edges, implying that the excitons have to undergo more 
relaxation steps to be trapped. Second, the bias-induced 
localization of the exciton states decreases their probabil- 
ity overlap *^2 n c 2 n c 2 , n , thereby reducing the relaxation 
rates W ss > . 

The above observations form the basis for understand- 
ing the results in Fig. [3J For small quenching amplitudes, 
r <C Wo, the downward relaxation is the fastest process 
and the trapping step becomes rate-limiting. Eventually, 
most of the population will end up in the exciton states 
near the lower band edge, where it waits to be trapped. 
In the absence of a bias (A = 0), these states have negli- 
gible quenching rates T s , as explained above, leading to 
large r values. After the bias is turned on, the amplitudes 
of the band edge states at the trapping site n = 1 steeply 
increase, giving rise to a strongly enhanced quenching 
rate. This will in turn result in a decrease of the overall 
harvesting time r. For very large biases, the overall har- 
vesting time will eventually rise again, as a result of the 
reduced relaxation rates W ss i . All this is in agreement 
with the behavior observed in Fig. [3J for V <C Wo- 

The nonmonotonic A-dependence of r observed in 
Fig. [3J (reaching a minimal t value at small bias) oc- 
curs when the quenching and scattering amplitudes are 
comparable, T ~ Wo- This nonmonotonic behavior is 
caused by the competition between the bias-induced ef- 
fects discussed above. More specifically, for small but in- 
creasing values of A, the increase in the quenching rates 
F s for the exciton states near the bottom of the band 
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FIG. 5: Harvesting time r [defined in Eq. (JSJ] versus the 
applied bias magnitude A calculated for a disorder-free chain 
of 50 sites. The solid lines correspond to the model where 
coherences are taken into account, while the dashed lines only 
account for populations and are identical to the results in 
Fig. E] 

is the dominant effect. This explains why r initially de- 
creases as a function of A. A further increase in A will 
shift the redistribution to the lower band edge, while the 
P s (0)-distribution shifts to the upper band edge, making 
the conditions for quenching less favorable. The connec- 
tion of these two distributions via the intra-band relax- 
ation accounts for the reduced harvesting efficiency be- 
yond some optimum bias value, due to the decrease of the 
relaxation rates for increasing A. Thus, we conclude that 
at r ~ Wq , there exists an optimal magnitude of the bias 
that provides the most efficient harvesting conditions. 

It should be noted that for large values of the bias, 
when the exciton states have become strongly localized, 
r becomes independent of the quenching amplitude, be- 
cause the harvesting process is limited in rate totally by 
the intraband relaxation. This is confirmed by Fig. [3l 

Figure [5] shows the corrections that a crude model (see 
Sec. IIII D[) for coherent transport yields, where the time 
evolution of the coherences is assumed to be simply expo- 
nential. The corrections are at most somewhere around 
25 percent, which makes sense as direct quenching only 
accounts for a fraction of the total quenched population. 
Also, dephasing of the coherences is provided through 
both the exciton-phonon coupling and the trap, and co- 
herent effects are therefore suppressed for all but the 
fastest energy transfer processes. A full discussion of the 
effects of coherent transport in unbiased chains should 
not be limited to the secular treatment provided here, 
and will be provided in a later paper. 

B. Disordered biased chain 

In this section, we present a study of the bias- 
dependent harvesting time r in disordered chains. We 




FIG. 6: Effect of uncorrelated on-site Gaussian disorder of 
strength a on the bias dependence of the harvesting time r 
calculated for two quenching amplitudes F = 10 _1 J (upper 
panel) and F = 10 J (lower panel). 



considered on-site Gaussian disorder (uncorrelated for 
different sites) with zero mean and standard deviation a, 
referred to as the disorder strength. Figure [6] shows the 
results obtained for two quenching amplitudes T — 1CP 1 J 
and T = 10J, keeping Wo = 10J and N = 50 in all sim- 
ulations. For other values of the scattering or quenching 
amplitude, similar results were obtained. The general 
observation that can be deduced from Fig. [5] is that the 
inclusion of disorder leads to a reduced efficiency of the 
energy harvesting process. The effect is especially sig- 
nificant for small A. Another consequence of disorder is 
that, in the case of T ~ Wo, the magnitude of the bias 
that is optimal for harvesting shifts to larger values. 

The explanation of these observations is straightfor- 
ward. As we already mentioned, the disorder results 
in Anderson localization of all the states on segments 
of the chain. First of all, this localization breaks the 
correlation that exists at low magnitudes of the bias 
(NA -C -Bjv - Ei = 4. 207 J) between the distributions 
P s (0) and r s : the states that are mostly excited initially 
and those that are quenched efficiently do not overlap 
anymore. Secondly, the localized states that have a sig- 
nificant amplitude at the trap are now not necessarily 
the lowest energy states and can hardly be reached dur- 
ing the downward intra-band relaxation. In addition, the 
localization slows down the overall exciton relaxation, 
because of blocking of the diffusion over the localized 
states at zero temperature.— All these factors result in 
a reduced harvesting efficiency at zero bias. Higher bias 
strengths tend to (partially) restore the correlation be- 
tween the energy and location of an exciton state, in 
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particular the lowest-energy exciton states will tend to 
lie in the vicinity of the trap. As a result, the states that 
are quenched most strongly are more accessible through 
intraband relaxation. This is why for small bias values 
A, the harvesting time r always decreases for increasing 
A, independently of the quenching amplitude T. 

For larger bias magnitudes, such that NA > En — 
Ei = 4.207J, the dynamic localization dominates over 
the disorder-induced localization. As a consequence, in 
this region the effect of introducing disorder is very lim- 
ited, in full agreement with the results presented in Fig. [5] 
Note that, like the Bloch localization discussed in the 
previous section and for similar reasons, Anderson local- 
ization expands the range of validity of our model. 

V. SUMMARY AND CONCLUDING REMARKS 

We studied the efficiency of harvesting of excitation 
energy in chains of strongly coupled dipolar units, in the 
presence of a linear energy bias. The harvesting com- 
bines the excitation transport following initial excitation 
at one end of the chain and quenching of the excitation 
by a trap at the other end. The energy gradient was set 
towards the trap. We considered the quenching rate of a 
particular exciton state to be proportional to its occupa- 
tion probability at the position of the trap. Within this 
model, the energy harvesting efficiency is governed by the 
interplay of the direct quenching of exciton states due to 
their overlap with the trap and the exciton relaxation 
to those states which are quenched most efficiently. We 
found a complicated scenario for the energy harvesting 
efficiency as a function of the bias magnitude A. Most 
importantly, we found that a bias does not necessarily 
increase this efficiency. 

The effect of a bias strongly depends on the ratio be- 
tween the amplitudes T (efficiency of the trap) and Wq 
(vibration- induced scattering). In the limit of T <C Wq, 
the initial population first scatters toward the bottom of 
the band, where quenching (the slowest process) eventu- 
ally takes place. Introducing a bias now greatly reduces 
r, because of its above mentioned tendency to shift the 
strongly quenched states to the lower band edge. Sup- 



pression of the exciton's downward relaxation upon in- 
creasing the bias magnitude now plays a less important 
role, because the relaxation is not the rate limiting step 
in the process. 

In the intermediate case of T ~ Wq, it depends on the 
bias strength what is the rate limiting step in the harvest- 
ing process. As a consequence, the bias dependence of 
the harvesting time r reveals a nonmonotonic behavior, 
in other words an optimal magnitude of the bias exists 
at which excitons are harvested most efficiently. As has 
been explained in Section [ill! our model is not expected 
to be appropriate in the limit of T ^ Wq , corresponding 
to what we expect is an unphysically strong quenching 
effect, as it does not adequately include coherent effects. 

We have also investigated how random energy disorder, 
superimposed on the overall bias, affects the above find- 
ings. We have found that disorder leads to a decrease of 
the harvesting efficiency, in particular for low bias magni- 
tudes, and gives rise to a nonmonotonic bias dependence 
of this efficiency, independently of the ratio T/Wq. We 
have shown that these effects can be understood from 
the disorder-induced (Anderson) localization of the exci- 
ton states. The effects of disorder are smeared at higher 
bias magnitudes, when the dynamic (Bloch) localization 
becomes dominant. 

To conclude, we note that the competition between the 
enhancement of transport and quenching, which under- 
lies the interesting nonmonotonic behavior discussed in 
this paper, is a very common theme in the physics of en- 
ergy harvesting systems. Although the model presented 
here is inspired by one-dimensional Frcnkcl exciton sys- 
tems, such as J-aggregates and conjugated polymers, it is 
to be expected that similar effects occur in other physical 
realizations, such as linear arrays of resonantly coupled 
quantum dots. 
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